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Abstract 



We investigate an algebraic model for the quantum oscillator based upon the Lie superalgebra 
5[)(2|2), known as the Heisenberg-Weyl superalgebra or "the algebra of supersymmetric quantum 
mechanics", and its Fock representation. The model offers some freedom in the choice of a 
i^H , position and a momentum operator, leading to a free model parameter 7. Using the technique 

C^ ' of Jacobi matrices, we determine the spectrum of the position operator, and show that its 

wavefunctions are related to Charlier polynomials C„ with parameter 7^. Some properties of 
these wavefunctions are discussed, as well as some other properties of the current oscillator 
model. 

> 

Q^ ■ 1 Introduction 

(N 

en ■ The relation between (one-dimensional) quantum oscillator models and orthogonal polynomials 

'sj" . (of a single variable) is very deep. For the canonical quantum oscillator, it can be found in any 
textbook that the position wavefunction is given in terms of Hermite polynomials. Among the first 
alternative oscillator model is the paraboson oscillator |251 I13j. where the position wavefunctions 
are given in terms of Laguerre polynomials. Depending on the oscillator model under consideration, 
the spectrum of the position operator can be continuous or discrete (with finite or infinite support), 
and so the wavefunction can correspond to a continuous or discrete orthogonal polynomial. Another 
C^ . example of an oscillator model with a continuous position spectrum is the su(l, 1) oscillator studied 
by Klimyk [20j, where the wavefunctions are given in terms of Meixner-PoUaczek polynomials, or 
its deformation the su{l, 1)^ oscillator [TB] with wavefunctions given in terms of continuous dual 
Hahn polynomials. Note that for all of the examples mentioned there is an algebra underlying 
the model, namely the Heisenberg-Weyl algebra, the paraboson algebra osp(l|2), the Lie algebra 
su(l, 1) and its deformation su(l, 1)^, respectively. 

In recent years there was also more interest in so-called finite or discrete oscillator models, 
inspired by applications in quantum optics or signal analysis. The seminal paper in this area is [3], 
where an oscillator model based on the Lie algebra 5u(2) was given, with position wavefunctions 
expressed in terms of symmetric Krawtchouk polynomials. This model was extended by an extra 
parameter in |141I15]. with an underlying deformed su(2) algebra, and position wavefunctions given 
by means of Hahn polynomials. 

Lately, also Lie superalgebras came into the picture as algebras underlying oscillator models. 
In [l7], the Lie superalgebra s[(2|l) and its finite-dimensional representations were used to give 
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an oscillator model with general Krawtchouk polynomials as wavefunctions (so with a discrete 
and finite position spectrum). Following this, infinite-dimensional discrete series representations of 
s[(2|l) were constructed in p[8], and this gave rise to an oscillator model with Meixner polynomials 
as wavefunctions (with a discrete and infinite position spectrum). 

In the current paper, we manage to construct an oscillator model for which the position wave- 
functions correspond to Charlier polynomials. The algebraic context is provided by the Heisenberg- 
Weyl superalgebra 5f)(2|2), an algebra generated by one boson and one fermion pair. This algebra 
is usually cited in the framework of supersymmetric quantum mechanics. 

In the following section, we shall introduce the superalgebra sf)(2|2), and more precisely a 
subalgebra S of its universal enveloping algebra. All (anti) commutation relations of this algebra 
are given, together with the action of its basis elements on the Fock space representation V of 
5f)(2|2). With the purpose of describing an oscillator model, we study the spectrum and the formal 
eigenvectors of a general self-adjoint odd element g of 5 in section 3. This is done in a broader 
context, since there has been some interest in the general construction of so-called superalgebra 
eigenstates [2]. The spectral analysis of this operator q leads to a Jacobi matrix, and yields the 
infinite discrete spectrum §={..., — -y/S, — \/2, —1,0, 1, \/2, ^/S, . . .} for q, with eigenvectors given 
in terms of Charlier polynomials. In section 4 we turn to the sf)(2|2) oscillator model: this is an 
algebraic model for which the Hamiltonian, the position and the momentum operators are elements 
of 5f)(2|2) (or S). The position operator is just the operator q studied in section 3. So we do 
indeed come up with an oscillator model for which the position wavefunctions are given by Charlier 
polynomials. In that section, we also present some properties of plots of these wavefunctions, 
and discuss a limit relation to another oscillator model. Since the momentum eigenvectors are 
also easy to give for the 5p)(2|2) model, it is a simple task to construct the corresponding sf)(2|2) 
Fourier transform kernel explicitly. The section ends with the computation of some other physical 
quantities of the model, such as uncertainty relations. In the final section, we give a differential 
operator realization of the superalgebra S, and end with some conclusions. 

It should be noted that Charlier polynomials have been related before to the Weyl algebra (and 
thus to the oscillator), e.g. in [ITIET]. In particular, they appear as overlap coefficients between 
eigenvectors of a number operator and a "shifted" number operator [27]. In the current paper they 
are, to our knowledge, for the first time interpreted as position wavefunctions of an oscillator. 

Finally, in this introduction we have spent some emphasis on the relation between algebraic 
oscillator models and related orthogonal polynomials appearing as wavefunctions. We should of 
course mention that also in the context of q-deformations of oscillator models and g-orthogonal 
polynomials, a large amount of work has been performed, see [19] and references therein. 



2 The Lie superalgebra S of the supersymmetric oscillator and its 
Fock representation 

The supersymmetric version of the one-dimensional harmonic oscillator has been studied by many 
authors, mainly in the context of supersymmetric quantum mechanics |281 [29l [9l \T0[ [5]. In [5l [6], 
the emphasis is on determining superalgebras of symmetries and supersymmetries, and relations to 
known Lie superalgebras have been studied. Let us reconsider here the case of a one-dimensional 
supersymmetric oscillator, the simplest example of A^ = 2 supersymmetric quantum mechanics. 
Such an oscillator is realized in terms of a bosonic creation and annihilation operator b^ with 
commutation relation 

[b-,b+]=l, (1) 

and a fermionic creation and annihilation operator a^ with anticommutation relations 

{a~,a'} = {a+,a'^} = 0, {a~,a+} = 1. (2) 



Furthermore, the two sets of operators commute with each other: 

[b^,a^] = 0, e,^ £{+,-}• (3) 

The Lie superalgebra generated by the even elements l,b^,b~ and the odd elements a'^,a~ is 
known as the Heisenberg-Weyl superalgebra sf)(2|2). The algebra of symmetries or the invariance 
superalgebra of the supersymmetric oscillator is usually identified as a subalgebra of the enveloping 
algebra W(sf)(2|2)). Let us consider here the Lie superalgebra S C W(sf)(2|2)) with four odd basis 
elements 

F+ = a+, F-=a-, Q+ = b+a-, Q- = b-a+, 

and even basis elements 

E+ = 6+, E- = b-, H = b+b- + a+a^- and L 

As a supersymmetric quantum system, H is the Hamiltonian and Q are the supercharges com- 
muting with H and satisfying {Q~^,Q^} = H. The complete set of (anti)commutation relations 
for the odd elements of S is given by: 

{F±,F±} = {Q±,Q±} = 0, {F+,F-} = 1, {g+,Q-} = if, 

{F±,Q^} = 0, {F+,Q+} = E+, {F',Q-] = E'. (4) 

The even-even commutation relations read 

[E-,E+] = l, [H,E^] = ±E^, (5) 

where of course 1 is a central element commuting with all elements. Finally, the even-odd commu- 
tation relations are 

[E^,F^] = [E^,F^] = 0, [E^,Q^] = 0, [E^,Q^] = tF±, 

[F,F±] = ±F±, [H,Q^] = 0. (6) 

The common Fock representation V of sf)(2|2), generated by a vacuum vector |0) and by the 
relations 

6-|0) =a"|0) = 

carries also an irreducible representation of the Lie superalgebra S. An orthonormal basis of V is 
given by the vectors 

^ ' -{a+y\0), me {0,1,2,...}, j e {0,1}. 



It will be convenient to write these vectors as \n), with n = 0, 1, 2, . . ., so y can be identified with 
the Hilbert space f{Z+): 

|2m) = i^|0), |2m + l) = i^a+|0). (7) 



Now it is straightforward to determine the action of the odd elements of S on this basis: 
F^\n) = \{l + (-inin + 1), F-\n) = 1(1 - (-iDIn - 1), 



Q^\n) = 1(1 - (-iDA/^ln + 1), Q'\n) = 1(1 + (-1)") Jf 1^ - 1). (8) 



Similarly, one finds for the action of the even elements: 



E-'W) = \[{l + (-iDi/^ + (1 - (-in V ^ ]\n + 2), 




This Fock representation is unitary with respect to the ^-structure 

(F±)t = F^, (Q±)t = QT^ (^±)t = ET^ ^t = H^ 



(9) 



(10) 



which follows from the common conjugacy relations for the boson/fermion creation and annihilation 
operators h^ and a^. 



3 Eigenstates of a self-adjoint odd element of S 

The algebra eigenstates associated to a real Lie algebra have been introduced by Brif [8]; they can be 
considered as generalizations of coherent or squeezed states. In [2], some superalgebra eigenstates 
of the Heisenberg-Weyl superalgebra sf)(2|2) have been computed, and properties of these states 
were investigated. In the current section, we shall determine the eigenstates (and eigenvalues) of a 
general self-adjoint odd element of the real Lie superalgebra S. Given the conjugations (fTOjl . the 
most general self-adjoint odd element, up to an overall factor, is 



7F+ + Q+ + 7F- + Q-, 



[ID 



where 7 is some real number. In the ordered basis {|n), n = 0, 1, 2, . . .} of the Fock space V , the 
operator q is represented by an infinite symmetric tridiagonal matrix: 
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(12) 



/ 



The form of this matrix follows from the actions dS])-®. For 7 > 0, such a matrix is a Jacobi 
matrix, and its spectral theory is related to orthogonal polynomials [TJ [23l I22j . It will soon be 
clear that only 7^ plays a role in the eigenvalues and eigenvectors of g, so let us assume for the 
moment that 7 > (the case 7 = is trivial, since the matrix (|12|) decomposes into irreducible 2x2 
blocks). The procedure to deal with such matrices is described in [23l §2]. One should construct 
polynomials Pn{x) of degree n in x, with p^i{x) = 0, Pq{x) = 1, and 



XP2n{x) = \/np2n-l{x) +-tP2n+l{x), 
Xp2n+l{x) = lP2n{x) + \/ U + lp2n+2{x) , 



(n = 0,1,2,...). 



(13) 



Such polynomials are orthogonal for some positive weight function w{x), and the spectrum of q is 
the support of this weight function. There is a mathematical condition related to this: the so-called 
(Hamburger) moment problem for the Jacobi matrix should be determinate [3 [22]. For a general 
tridiagonal Jacobi matrix, a sufficient condition [TJ [22] for this is that the (infinite) sum of the 
matrix elements above the diagonal should be infinite, which is clearly satisfied for ()12p . So the 
spectrum of q is just the support of the weight function w{x). Furthermore, for a real value x 
belonging to this support, the corresponding formal eigenvector of q is given by 

oo 

v{x) = '^Pn{x)\n). (14) 

n=0 

The recurrence relations (I13p lead to a solution of the polynomials Pn{x) in terms of terminating 
hyper geometric series; for their notation we follow that of standard books [Il[31[26]. The polynomials 
that will appear here are the Charlier polynomials Cn{x]a), defined by [211 [121 [1]: 

Cn{x;a)= ^fJ~'''J''--^. (15) 

These polynomials satisfy a discrete orthogonality relation: 

oo ^ 

V —rCm{x; a)Cn{x; a) = a-^'e^'nldmn (16) 

^-^ x\ 

when a > 0. Using the known forward and backward shift operator relations [21] 

C„(x; a) - Cn{x - 1; a) = Cn-i{x - 1; a), 

a 

X 

Cnix; a) Cn{x - 1, a) = C„+i(x; a), (17) 

a 

it is easy to see that the following holds: 

Proposition 1 For 7^0, the solution of the recurrence relations (|13|) is given by 

(—7)" 
P2n{x) = — ^=^C„(x^;7^), 



P2n+i{x) = - ^ ^A xCn(x^-l;7')- (18) 



Using the orthogonality relation of Charlier polynomials (J16p . one can deduce the corresponding 
orthogonality relation for the polynomials Pn{x)- For Charlier polynomials, the support is discrete 
and running from to +00; but due to the relations (jlSp the support for the p„'s is from —00 to 
+00 (also discrete). 

Proposition 2 The polynomials Pn{x) satisfy a discrete orthogonality relation: 

^^w{x)pn{x)pm{x) = e^ 6mn, (19) 

where 

S = {±Vk I keZ+} = {..., -\/3, -\/2, -1,0, 1, \/2, \/3, . . .}, (20) 

and where the weight function is given by 

<^) = Yw: fo^^ = ^^ (A; = 1,2,3,...) (21) 

and by w{x) = 1 for x = 0. 



So the spectrum of the operator q is discrete and given by (j20p . 

For further use, it will be convenient to introduce the corresponding orthonormal functions 
defined by Pn{x) = e~^ ''^■s/w{x)pn{x), i.e. 



P2n(x) = (-!)"■ -^^ ^^e ^ '^ ; :Cn(x^;T) 



^.2_ 



P2n+i{x) = S=-e-^ '^-L^==xCn{x^ - l'^^), (22) 



V2 v^n!(x2)! 

satisfying 



^Pm(a;)Pn(a;) = 5mn- (23) 

Following (|14p . the normalized eigenvectors of g, for an eigenvalue x S S, are given by 

00 
){x) = Y,Pn{x)\n). (24) 



00 
v( ] 

n=0 



4 A discrete oscillator model related to 5f)(2|2) 

In recent years, there has been some interest in the development of new oscillator models, partic- 
ularly in the context of finite oscillators. For these models, one requires the same dynamics as for 
the classical or quantum oscillator, but the operators corresponding to position, momentum and 
Hamiltonian can be elements of some algebra different from the traditional Heisenberg (or oscil- 
lator) Lie algebra. In the one-dimensional case, there are three (essentially self-adjoint) operators 
involved: the position operator q, its corresponding momentum operator p and the Hamiltonian H 
which is the generator of time evolution. The main requirement is that these operators should sat- 
isfy the Hamilton-Lie equations (or the compatibility of Hamilton's equations with the Heisenberg 
equations): 

[H,q] = -ip, [H,p]=iq, (25) 

in units with mass and frequency both equal to 1, and h = 1. Contrary to the canonical case, the 
commutator [q,p] = i is not required. Apart from ()25p and the self-adjointness, it is then common 
to require the following conditions [3]: 

• all operators q, p, H belong to some Lie algebra or Lie superalgebra A; 

• the spectrum of H in (unitary) representations of A is equidistant. 

It is within this framework that the oscillator models mentioned in the Introduction fit. In the 
current section, we shall examine the oscillator model related to s[](2|2) and describe some of its 
properties. 

4.1 The s{)(2|2) oscillator model and its position and momentum wavefunctions 

Clearly, it is possible to construct a new oscillator model using the superalgebra sf)(2|2) (or, more 
precisely, S) and its representation V . It might look a bit strange to build an ordinary oscil- 
lator model from a supersymmetric model, but at least it is possible to satisfy all the previous 
requirements within this model. 



For the Hamiltonian H^ one can in fact choose a diagonal operator with the same spectrum 
as the canonical oscillator. For this purpose, one extends the superalgebra sf)(2|2) by a parity or 
reflection operator fi, with action 

R\n) = (-l)"|n). 



Taking then as Hamiltonian 
one finds indeed the spectrum 



H = 2H + -R (26) 

H\n) = (n -\ — )\n). 

As in |18j . the position operator q is (up to an overall factor) the most general self-adjoint odd 
element of S. In other words, this is the operator (jlip q = jF'^ + Q^ + ^F~ + Q~ investigated in 
the previous section. Then the momentum operator must take the form 

p = i7F+ + iQ+ - i^F- - \Q- , (27) 

and ([2^ is satisfied. 

In the previous section, we have determined the spectrum S of the position operator q and its 
formal eigenvectors v{x) for x G S. The Jacobi matrix corresponding to p is very similar, and one 
obtains the following result: 

Proposition 3 The spectrum of the momentum operator p is also given by S. For an eigenvalue 
y G S, the normalized eigenvector u{y) is given by 

oo 

u{y) = Y,i^Pn{y)\n), (28) 

n=0 

where the expressions pn are the same as before, given by (I22p . 



We can now continue to investigate some properties of the sf)(2|2) oscillator model. First of 
all, the position (resp. momentum) wavefunctions of this oscillator are the overlaps between the 
(normalized) g-eigenvectors (I24p (resp. p-eigenvectors (I28p ) and the if -eigenvectors. Clearly, it is 
sufficient to consider only the position wavefunctions, denoted by Lpn {x) and given by 

^^;^\x)=pn{x), (29) 

where x belongs to S. The dependence of the wavefunctions on the model parameter 7 is made 
apparent in the notation. So we are dealing with a discrete position wavefunction (but with infinite 
support). Wavefunctions of this type have appeared also for the (infinite discrete) s[(2|l) oscillator 
model |18j . In the current case, the shape of these discrete position wavefunctions is very similar. 
We have plotted some examples of ipn {x) in Figure 1, for certain values of 7. Note that the shapes 
are indeed discrete versions of the common oscillator wavefunctions in terms of normalized Hermite 
polynomials. 

As far as the plots in Figure 1 are concerned, keep in mind that we are dealing with discrete 
wavefunctions with infinite support S. So only the dots in Figure 1 are part of the wavefunctions, 
not the curves connecting these dots (these curves have been drawn just to emphasize the "shape" 
of the discrete wavefunctions). Note also that we have plotted only a finite part of the infinite 
support (namely only the dots corresponding to the support [—vk,+vk\ with k = 10). 



4.2 A limit related to the position wavefunction 

It is not straightforward to consider direct limits of the 5f)(2|2) position wavefunctions ifn {x), 
e.g. when the model parameter 7 varies. Instead, we need to look for another oscillator model of 
which the current model itself is a limit. This is provided by the finite oscillator model related 
to the Lie superalgebra 5[(2|1) |17j . Recall that this s[(2|l) oscillator model has two parameters: 
a representation parameter j (which is a nonnegative integer) and a model parameter p (with 
< p < 1). The representations in [T7] are finite-dimensional (with dimension 2j + 1). The 
support of the s[(2|l) position wavefunctions is given by ztyfe (where A; = 0, 1, . . . , j). The actual 
5l(2|l) wavefunctions are given by [El (44)-(48)]: 



X = ibVK, A; = 0, 1, ... j; 



<^^r(x) = (-ir^/^^i^„(x2;p,j), (n = 0,l,...,j), 



X = ibVK, /c = 0, 1, . . . j. 



.2^+1 (x) = ±{-lY^^^Kn{x' - l;p,j - 1), (n = 0, 1, . . . , j - 1), 



Herein, Kn{x;p,N) is the normalized Krawtchouk polynomial |171 (24)], i.e. the Krawtchouk poly- 
nomial Kn{x;p,N) [21, 9.11] multiplied by the square root of its weight function and the inverse 
of its norm. 

In general, the representation parameter j and the parameter p are independent in the s[(2|l) 
oscillator model. Let us, however, consider the following dependence: 

p = -f^/j, (30) 

where 7 is a nonzero constant. Then the classical limit from Krawtchouk polynomials to Charlier 
polynomials [21, 9.14] implies: 

lim<f>^f/^'^\x)=^^;y\x). (31) 

In other words, for large values of j the position wavefunctions of the s[(2|l) finite oscillator model 
tends to those of the sfi(2|2) oscillator model, provided the parameter p is represented as 7^/j in 
the limit. 

The behaviour of these wavefunctions for large values of j is also confirmed in plots. In Figure 2, 
we have plotted the ground state (n = 0) and first excited state (n = 1) of the 5[(2|1) oscillator for 
increasing values of j, namely j = 5, 10, 20, 30 (note that from j = 20 onwards, not the complete 
finite wavefunction is plotted, since we plot the support only in the interval [—4,4]). For these 
plots, we have taken 7 = 1 as an example. In the bottom line of Figure 2, we have plotted the 
corresponding wavefunctions of the sf)(2|2) oscillator, again for the value 7=1 (and again in the 
interval [—4,4]). From these plots, the limit relation is clearly illustrated. 

Recall from |17] that for large j-values the behaviour of the s[(2[l) wavefunctions 4>n' (x) is 
the same as the behaviour of paraboson wavefunctions ^n " (x) for large values of a. In the 
process just described, however, p is a fixed constant. This explains why it is not possible to give a 
direct limit from the s[(2[l) wavefunctions to a known family of continuous wavefunctions like the 
paraboson wavefunctions or the canonical oscillator wavefunctions. 



4.3 An expression of the s[)(2|2) Fourier transform 

In canonical quantum mechanics, the momentum wavefunction (in L^(M)) is given by the Fourier 
transform of the position wavefunction (and vice versa), with kernel K{x,y): 



i^iy) = I K{x,y)4>{x)dx, K{x,y) 



1 



^-ixy 



This means that K{x, y) is also the overlap of the (formal) position eigenvector for the eigenvalue 
x with the momentum eigenvector for the eigenvalue y. In the current case, the kernel for the 
corresponding sf) (2 1 2) Fourier transform is 



K^''\x,y) = {Hx),uiy)) = X;(-i)>(7)(x)^W(y), 



(32) 



n=0 



where ipn (x) = Pn{x) is given by (i22|) . 

It is easy to compute this kernel explicitly. Splitting the sum over even and odd n- values, 
using (I22p . and using the known bilinear generating function for Charlier polynomials |24j . 



— r 

V — C„(x; a)Cm{y; a) = e\l - t/af+y ^Fo 



m=0 



-x,-y 



t 



(a-i)V' 



(33) 



one finds 



K(^\x,y) 



-'^i'^ lo^W+y'^ 



■{2^y 



2V(x2)!(y2)! 

xy 



iY 



2-Pb 



l-x2, 1-2/2 



1 



-x2,-y2 



1 

'4^ 



472 



(34) 



Herein, x and y are elements of S, i.e. they are of the form x = :izVk, y = ±Vi with k and / 
nonnegative integers. 

4.4 Some other physical quantities of the sf)(2|2) oscillator 

As we are speaking of an oscillator model, let us compute here some other quantities that play a 
role in physics (all under the assumption m = u = H = 1). Taking into account the expression for 
the position operator q (jlip and the momentum operator p (127p . one finds 



[g,p] = -2i(7'[F+,F-] + [Q+,Q-]). 
Using the actions ([8]) in the Fock space, this leads to (n = 0, 1, 2, . . .) 



(-ir 7 



-n + 
2 ' 



1 -(-!)" 



[q,p] \n) = 2i 

or for even and odd states separately: 

[q,p]\2m) = 2i{'y — 112)12171), 
[q,p]\2m - 1) = -21(7^ - m)\2m - 1). 



n> 



(35) 



These relations are similar to the ones for the paraboson oscillator [13, (16)], apart from the m- 
dependence. 



Another quantity that can easily be computed is the standard deviation for q and p when the 
system is in the stationary state |n). Using 

{Aq)n = \/{n\q'^\n) — {n\q\n)'^ 

and the expHcit action of g = 7F+ + Q^ + jF^ + Q^ , one finds 



The computation for (A]5)„ is the same, so one finds as "uncertainty relation" quantity 

(Ag)„(Ap)„ = 7^ + in + i(l-(-in. 

Although this depends on the state number, note that its value is the same for two consecutive 
states: 

"2 -2 

Also the operator ^ + ^ is worth considering. On the algebraic level, one finds 

following from the explicit expressions of p, q and dH). For the action in the Fock space, this implies 

(y + \)\n) = (7' + 2 + ^^)l^)- 

Observe that all these quantities yield minor deviations from the corresponding ones for the 
canonical oscillator or for the paraboson oscillator. 

5 Some concluding remarks 

Before coming to some conclusions, let us mention a differential operator realization of the operators 
studied here. It is well known that the Heisenberg-Weyl superalgebra and its Fock representation 
can be realized in terms of an ordinary variable z and a Grassmann variable 9 (with 9^ = 0). The 
Fock representation states are then given by 

z" z" 

|2n) = ^=, \2n + l) = ^=9. (36) 

\/n\ yn\ 



The Lie superalgebra basis elements take the form 

F+ = 9, F- = de, Q+ = zde, Q- = 98,, 

E+ = z, E- = d,, H = zd, + 9de. (37) 

In this realization, the position operator corresponds to the differential operator 

q = 9{^ + d,) + de{i + z). (38) 

It is worthwhile to give the eigenfunctions of this operator explicitly. This can be done using the 
normalized eigenfunctions ()24p of q and by using the generating function for Charlier polynomi- 
als [a] 

Y^Cn{x;a)- = e\l--Y. 
10 



One finds: 

1 e-7V2-72 , , ^ 

v{x) = -^ ^—(7 + zf'Hl + z ± Vk9) 

forx = ±\/fc, A; = 1,2,3,... (39) 

and v{0) = e~'^ /2-7^_ j^i ^j^jg realization, it is easy to verify that qv{x) = xv{x), for x G S. 

To summarize, we have worked out a new model for a quantum oscillator based on the Lie 
superalgebra st)(2|2) and its Fock representation. The Hamiltonian H, the position q and the 
momentum p of the model are three self-adjoint elements of the enveloping algebra of sf)(2|2) 
satisfying the Hamilton-Lie equations (I25p . The model allows a free parameter 7 which appears in 
the expression of q. The spectrum of the energy operator H coincides with that of the canonical 
quantum oscillator. The spectrum of the position operator q is infinite discrete and given by § = 
{. . . , — v3, — v2, —1, 0, 1, v2, v3, . . .}. The most interesting fact is that the position wavefunctions 
^n (x) correspond to "symmetrized" Charlier polynomials (i.e. even wavefunctions are Charlier 
polynomials in x^, and odd wavefunctions are x times Charlier polynomials in x^). Some plots of 
the position wavefunctions ipn (x) reveal interesting properties. In particular, the wavefunctions 
can be seen as limits of finite discrete wavefunctions appearing in an s[(2|l) oscillator model. 

The current model is sufficiently simple, and could perhaps be extended to more dimensions, 
e.g. using the supersymmetric Heisenberg-Weyl algebra sf)(2|(i). Whether such an extension gives 
rise to interesting results related to orthogonal polynomials, is a question that will be studied in 
the future. 
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Figure 1: Plots of (part of) the wavefunctions (pn (x) for n = 0, 1,2,3 and for some values of 7, 
where x £ S. 
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Figure 2: Plots of the s[(2|l) wavefunctions (pn' (x) for n = 0, 1 with p = 7^/j and for some values 
of j (with 7 = 1), and their comparison with ipn (x) (bottom figure). 
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